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What is claimed is: 

1. A method for constructing a composite response surface based on neural 
networks and selected functions, the method comprising: 

(1) providing a set of h initial parameters that determine variation of 
provided data for a target variable, where each parameter corresponds to a 
coordinate in an h-dimensional parameter space G; 

(2) providing a decomposition of the h parameters into a first set of s simple 
parameters f s , numbered i = 1 , . . . , s, that may be used to describe the provided data 
with polynomials of total degree no greater than a selected number M s , and a 
second set of c complex parameters g j? numbered j = 1, ... , c, that may be used to 
describe the provided data using neural networks, and with s + c = h, where s, c 
and M s are selected positive integers; 

(3) providing a simplex, having s+1 vertices, numbered k = 1, ... , s+1, and 
centered at a selected point in the space G; 

(4) applying a neural network for each of the s+1 vertices, and training each 
of the s+1 neural networks, using selected simulation data obtained by varying the 
parameters gj to generate a first sequence of network functions R k (g l5 ... , g c ); 

(5) providing a second sequence of shape functions P k (f 1? ... , f s ) that satisfy 
the conditions P k (fi, . . . , f s ) = 1 at the vertex numbered k and P k (f„ . . . , f s ) = 0 at 
any vertex other than vertex number k, and £ P k (f l5 . . . , f s ) = 1 for all values of f„ 
... , f s ; and 

(6) forming a composite function CRS(fj gj) defined by 

s+1 

CRStf, gj } = S P k (f„ ... , f s )R k (g„ ... , g c ). 
k=l 
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2. The method of claim 1, further comprising selecting said set of complex 
parameters to include at least one polynomial in said complex parameters gj having 
a selected degree M c satisfying M c > M s . 

3. The method of claim 1, further comprising choosing said integer M s from 
the group of integers consisting of 1, 2 and 3. 

4. The method of claim 1, further comprising selecting said set of complex 
parameters to include any of said h parameters that does not qualify as a simple > 
parameter. 

5. The method of claim 1, further comprising: 

(7) providing an objective function OBJ(f k ,gj) n , dependent upon at least one 
of the parameter values f v . . . , f s , g v . . . , g c , for the composite function CRS{f k ,gj} 
at each of N selected locations in G space, numbered n = 1, . . . , N, associated with 
the target variable, and providing a corresponding objective function value OBJ n 
for the target variable at each of the N selected locations, where n is a selected 
positive integer; 

(8) computing a training error value TE{gj,} as a non-negative weighted sum 
of functions of differences F n (OBJ n - OB J(f k ,gj) n ), where each function F n is 
monotonically increasing in a magnitude of the function argument and has a value 
0 where the function argument is 0; 

(9) when the training error value TE{gj} is greater than a selected threshold 
error value e, providing at least one of a modified set of shape functions P k (f l5 

f s ), and returning to step (6); and 
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(10) when the training error TE{gj} is no greater than the threshold error 
value e, accepting the present composite response surface. 

6. The method of claim 1, further comprising applying said composite 
response surface to optimization of a design of a physical object. 

7. The method of claim 6, further comprising choosing said physical object 
to be a shape for an aircraft component. 

8. The method of claim 1, further comprising applying said composite 
response surface to modeling of a response to a process. 

9. The method of claim 1, further comprising applying said composite 
response surface to modeling response of a physical object. 



